Attosecond Transient Absorption Spectroscopy of doubly-excited states in helium 
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Strong-field manipulation of autoionizing states is a crucial aspect of electronic quantum control. 
Recent measurements of the attosecond transient absorption spectrum of helium dressed by a few- 
cycle visible pulse [Ott et al., arXiv: 1205.0519 physics. atom-ph] ] provide evidence of novel ultrafast 
resonant phenomena, namely, two-photon Rabi oscillations between doubly-excited states and the 
inversion of Fano profiles. Here we present the results of accurate ab-initio calculations that agree 
with these observations and in addition predict that (i) inversion of Fano profiles is actually periodic 
in the coupling laser intensity and (ii) the supposedly dark 2p^ state also appears in the spectrum. 
Closer inspection of the experimental data confirms the latter prediction. 

PACS numbers: 31.15.ac, 32.80.Fb, 32.80.Rm, 32.80.Zb 
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The fast development of attosecond laser pulses in 
the last decade gave access to the time-resolved study 
of correlated electron dynamics in atoms, molecules and 
solids on their natural time scale In this context, 
attosecond transient absorption spectroscopy (ATAS) [H 
is emerging as a prominent technique complementary to 
the technologies which detect charged photofragments, 
like COLTRIMS Q and VMI H, to monitor and control 
transiently bound states. Compared to photofragment 
detection, ATAS provides higher energy resolution, and 
is applicable to condensed phases Q as well as to the gas 
phase, thus making it a good candidate for the investiga- 
tion of ultrafast electron dynamics in chemically relevant 
samples. Moreover, on the theoretical side, ATA spectra 
are significantly less demanding to reproduce than pho- 
toelectron distributions. While the latter require specifi- 
cation of scattering boundary conditions, the former only 
need the electronic wave function in the vicinity of the 
reaction center. 

So far, ATAS has been employed to track both field- 
free ultrafast dynamics of coherent ensembles of valence- 
hole excited states 0, 0] , as well as the driven dynam- 
ics induced by strong external visible (VIS) or near- 
infrared (NIR) dressing fields 04l3 , including the follow- 
ing processes: electromagnetically induced transparency 
(EIT) |i,[l^; instantaneous ac-Stark shift 0] ; multipho- 
ton transition matrix elements [l2j ; interference of wave- 
packets [ll|; and time- resolved autoionization 0. 

In a recent experiment [§], ATAS has been used to 
investigate helium under the influence of a VIS dress- 
ing field at tunable intensity in the energy region of the 
doubly-excited states converging to the N=2 threshold. 
The experiment revealed clear Autler-Townes splitting 
of the 2s2p autoionizing state (the lowest term in the 
SP22 bright series) due to the coupling between this 
state, the ^S state, and the consecutive spjg state, 
and inversion of the characteristic Fano profile asymme- 



try in the higher terms of the spj„ series as a function 
of both the pump-probe time delay and of the intensity 
of the VIS dressing field. 
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FIG. 1. Energy level scheme. The attosecond transient 
absorption spectrum of helium records the linear response of 
the atom, dressed by a strong external VIS field, in the ex- 
treme ultraviolet (XUV) energy range. A sub-femtosecond 
broadband XUV pulse excites a wide range of states in the 
single-ionization region. The population and phase of the 
states so created is altered as a result of the interaction with 
the strong dressing field, due to multiphoton or even non- 
perturbative transitions. The XUV dipole response of the 
atom refiects such composition changes. 



Here we present the results of accurate ab-initio calcu- 
lations that reproduce the experimental conditions in [§] . 
We compute the TAS of a sub-femtosecond extreme ul- 
traviolet (XUV) pulse for a helium atom dressed by a 
short (7fs) visible laser pulse (730 nm), as a function of 
the time delay between the two pulses and of the peak 
intensity of the dressing field. Our calculations are in ex- 
cellent agreement with most of the current experimental 
results and point towards further interesting quantum 
phenomena of the physics of two interacting electrons 
exposed to both the Coulomb field of the nucleus and a 
laser field. In particular, the inversion of Fano profiles 
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is actually found to be periodic in the laser intensity, in 
agreement with the interpretation of this effect in terms 
of the energy shift of the doubly-excited states, which is 
approximately proportional to the intensity of the dress- 
ing laser @. Moreover, clear signatures of the supposedly 
dark 2p^ state are observed. Finally, we show that 
the ATA spectrum can be construed in terms of the non- 
diagonal component of the linear electrical susceptibility 
of the dressed atom. 

In the ATAS experiment under examination, the weak 
60 eV XUV pulse £{t) induces a time-dependent dipole 
moment d{t;£) in the atom dressed by the VIS field, 
which in turn results in a dipole emission that coherently 
contributes to the XUV field, thus altering its spectrum 
(see Fig. [Ij. For an optically thin medium, the relative 
change of the spectral intensity of the XUV pulse can 
be expressed in terms of an effective transient absorption 
atomic cross section a{Lo;£) 
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where d and £ are the Fourier transform of the com- 
ponent of the dipole moment and the XUV field along 
the laser polarization, respectively, and where we em- 
phasize the functional dependence of the response on the 
full spectrum of the impinging attosecond pulse. In this 
process, strongly correlated doubly-excited states (DES) 
have a preponderant role. For the latter, common mod- 
els to evaluate the rate of non^erturbative processes, 
like ionization by tunneling [l3.[l5l|. are inapplicable, due 
to the complete breakdown of the single-active-electron 
(SAE), the adiabatic, and the non-polarizable parent- 
ion approximations 1^ 121- In this context, calculations 
based on simplified models d, Hsl-Ii^l, S AE, or station- 
ary methods, like Floquet theory [13|, |25[ have therefore 
limited applicability. The interpretation of these exper- 
iments require instead both a complete ab-initio repre- 
sentation of the system and a direct solution of the time- 
dependent Schrodinger equation (TDSE). 

Here, we solve the TDSE with a Krylov representation 
of a second-order exponential unitary time-step propa- 
gator in velocity gauge Reflection from the box 
boundaries is prevented by complex absorbing poten- 
tials. The wave function is expanded on a two-particle 
spherical basis, the angular part is represented by bipolar 
spherical harmonics [27[ and the radial part by B-splincs 
with an asymptotic spacing of 0.5 au. Each total an- 
gular momentum comprises all the partial-wave channels 
with configurations of the form N£ei' with N < 2, and 
a fuU-CI localized channel nin'i' that reproduces short- 
range correlations between the two electrons [2^. Such 
representations of the correlated electron basis is able to 
represent in a concise way and with high accuracy all the 
resonant series in helium converging to any given thresh- 
old 12811 .and, thanks to the localized channel, the first 



few doubly-excited states converging to thresholds not 
included in the close-coupling expansion as well. The 
ATAS response is evaluated with a velocity- gauge ana- 
logue of (dl) . obtained via an application of the Ehrcnfcst 
theorem 
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where A and p are the FT of the vector potential and 
the canonical momentum. The momentum p is given 
by the sum of two complementary terms: p-, computed 
numerically while in the presence of the dressing field, 
and computed analytically for the subsequent field- 
free evolution. The result is thus equivalent to the one 
that would be obtained for a simulation protracted for 
an infinite time. In the simulation, we include the states 
with orbital angular momentum up to £rnax = 5 for the 
localized channel, and a total angular momentum up to 
Lmax = 10. We checked for convergence with respect 
to both these limits by repeating few calculations with 
either L„iax ~ 15 or tmax = 10, as well as by including 
the N = 3 partial- wave channels, without observing any 
significant change in the result. 
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FIG. 2. ATAS as a function of tlie time delay, a) Theoretical 
transient absorption spectrum as a function of the XUV pho- 
ton energy (x-axis) and of the time delay (y-axis) between the 
XUV and the VIS pulses, for a VIS intensity 1=3.5 TW/cm^ 
b) Matching experimental spectrum recorded for a nominal 
VIS intensity of 2.5±1.5 TW/cm^. c) Theoretical, d) Exper- 
imental Fourier transform with respect to the time delay of 
the spectra shown in panel a-b). 

In Fig. [5]-a,b we compare the present theoretical ATAS 
spectrum (7(0;; r), as a function of the photon energy 
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and the time delay r, to the experimental optical den- 
sity (OD) (OD cx a). The agreement is excel- 
lent. The vertical features visible in the two spectra are 
due to the bright spj„ series of doubly-excited states, 
71 = 2, 3, . . . [29|. At large negative time delays, when 
the dressing pulse comes well before the XUV pulse, 
the ATAS spectrum displays the characteristic Fano pro- 
files 2^3lj of the field- free atom. When the two pulses 
overlap, the SP22 state separates in two branches; this 
is the well known 0, [l^ Autler-Townes splitting 
(ATS) associated with the strong coupling between the 
SP22 a-nd the 2p^ state, which, with the current en- 
ergy of the dressing field, 1.7 eV, is resonantly coupled 
to both the first and the second term of the spj„ se- 
ries (sec Fig. [Ij. The faint transversal wavelike features 
visible around both the SP22 the sp^g terms for pos- 
itive time delays are slices of the hyperbolic branches 
Tn{E) = 2TTn/\E — Eq\. These features are caused by 
the beating between the amplitude of the resonantly- 
enhanced two-photon transition from one XUV-excited 
autoionizing state to the continuum-coupled other res- 
onant state and the amplitude of the XUV transition 
from the ground state to the same coupled-resonant and 
continuum states. As far as the direct-ionization part is 
concerned, such phenomenon is the transient- absorption 
analogue to the one described for the photoelectron spec- 
trum in [32}. At very large time delays, the signal con- 
verges again to the field-free case. Here, the fringes con- 
verging around each resonance are a consequence of the 
sharp change in the phase or population of the local- 
ized part of the resonances induced by the strong VIS 
pulse. Finally, in the dark narrow region between the 
two Autler-Townes branches of the sp2 2 state the effec- 
tive transient absorption cross section is negative, what 
indicates that the XUV field in that region is actually 
amplified rather than absorbed, as a result of the inter- 
action with the atom. This observation confirms previous 
predictions [l3, 23, 33 1 and experimental findings Q- 

In Fig. [21-c,d we compare the FT with respect to the 
time delay of the ATA spectra in Fig. [2]-a,b [3J| 



[CT(a;,T) - cr(w,oo)] , (3) 



where Wr is the Fourier Energy. To interpret such bi- 
dimensional spectra, let us consider again Eq. [TJ Since 
the XUV pulse is weak, the dipole response of the laser- 
driven system can be assumed to be linear 



(icj'x(w, Lij')E{ui'). 



(4) 



where x(w, w') is the electric susceptibility of the dressed 
atom. For time-invariant samples, of course, the sus- 
ceptibility is diagonal: x(a;,a;') = 5{ijJ — uj')xd{u!). This 
is not the case here, however. Under the action of any 
given excitation frequency, a dressed atom will in general 
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FIG. 3. Detail of the comparison between theoretical (left 
panel) and experimental (right panel) ATAS as a function 
of the XUV photon energy (x-axis) and of the time delay 
between VIS and XUV pulses (y-axis) in the energy region 
between the spjj s-nd the spjg doubly-excited states. In the 
insets, closer inspection unveils a clear signature of the 2p^ 
resonance, which mediates the transition between the lowest 
bright states. 



give rise to a response at other frequencies as well. This 
means that, alongside the field-free diagonal component, 
the electrical susceptibility also comprises a non-diagonal 
component Xnd, 



X(w, w') = 5{uj - uj')xd{i^) + Xnd(w,w')- 



(5) 



For moderate intensities of the dressing lasers, the 
strongest features of Xnd are expected to be located close 
to the diagonals cj = w' ± 2nujyis. If the spectrum of the 
exciting laser is narrower than 2cj„i 
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as in \li 

transient absorption spectroscopy only probes the re- 
sponse close to the central diagonal u = uj' . This provides 
valuable yet incomplete information on the ultrafast opti- 
cal response of the dressed system. In the present ATAS, 
on the other hand, the spectral width of attosecond pulse 
exceeds 2ujms-, therefore the response at a; = iL:'±2ujyis re- 
sults in a detectable heterodyne signal. This result holds 
for isolated as well as for trains [lol - [l2l | of attosecond 
pulses. The variation of the XUV field due to the change 
in time delay r with respect to the dressing field, which 
defines the time origin, is readily parametrized as 



(6) 



where £0 corresponds to a pulse centered at the time 
origin. From Eqns. ([T][B]) and in the limit of extremely 
short single pulses, where fo(w) is eventually constant in 
an interval of the order of the excitation frequencies of 
the system, we obtain 



ic a{uj, oJt) 
(27r)3/2a; 



Xnd{uJ, UJ-UJr) - Xndi^^ + ^r) (7) 



The strong signals close to the ojr = axis, therefore, 
are associated to the electromagnetically induced trans- 
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parency (EIT) of the DES that is visible also in femtosec- 
ond TAS [13, 0, 24 [ . The two strong features above the 
SP22 and the spjg resonances at ojr ~ "^i^vis ~ 3.4 eV, 
instead, are heterodyne signals that are visible only in 
ATAS. In particular, the one above the spjg state can 
be assigned to the Xndi^sp+ ~ '^^^ 



) compo- 



nent of a, while the one above the spjg state to the 
2w^is,a;„„+ ) component. Both signals are 

strongly enhanced by the presence of the resonant 2p^ ^ S 



intermediate state (compared to the higher terms in the 
sp2n series) and comprise a non-resonant as well as a 
resonant part. The latter is indicative of an incipient res- 
onant two-photon Rabi oscillation between the two 
states (see Fig. [T]). 

As a consequence of its even parity, the 2p^ is nor- 
mally a "dark" state. Yet, due to the effect of the dress- 
ing field both on it and on the ground state, the dipole 
amplitude between these two states is not exactly zero. 
Indeed, two faint signals at photon energy ~ 61.5 eV 
and 62.3 eV close to the = axis are visible in 
Fig. [2] To confirm this observation, in figure [3] we show 
the theoretical and experimental ATAS with enhanced- 
contrast insets close to where the 2p^ is expected to 
emerge. The upper Autler-Towncs branch of this state 
is clearly visible in both insets with an intensity compa- 
rable to that of the narrow sp2^3 resonance. This 
result shows that TAS is not restricted to the study of 
bright states only. Also "dark" states which are coupled 
to bright ones can emerge. 

Fig. 2^ shows the experimental ATAS for overlap- 
ping XUV and VIS pulses, as a function of the in- 
tensity of the dressing laser. The spectrum is aver- 
aged over a time-delay window of 1.2 fs to reduce the 
noise. Fig.|3|DC show the theoretical prediction for a fixed 
time delay r = 200 au (maximum ATS). Overall, ex- 
cept for the highest intensities, the agreement is good: 
we observe the increase of the ATS of the SP22 reso- 
nance as well as the inversion of the asymmetry of the 
Fano profile for all the higher terms in the spj„ series 
as the intensity of the dressing field increases. In 
this latter effect has been attributed to the ac-Stark 
shift AEac-stark{t) of the localized component autoion- 
izing states. The energy shift translates into an extra 
phase A(/)(t) = ~ AEac-stark{t')dt' between reso- 
nant and background amplitude that is approximately 
proportional to the peak intensity of the laser. As a con- 
sequence, the asymmetry (Fano profile) inversion is ex- 
pected to be periodic in the laser intensity. Our calcula- 
tions (see Fig. 2]) clearly confirm this behavior. In the ex- 
perimental spectrum, all resonant features suddenly fade 
for intensities larger than 4 TW/cm^, close to the instru- 
mental limit, which is sufficient to observe just the first 
inversion for the highest terms in the series. Convolution 
of the theoretical results in Fig. \^ with the experimen- 
tal energy and intensity resolutions (cte ~20 mcV and 
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FIG. 4. ATAS for overlapping XUV and VIS pulses as a func- 
tion of the XUV photon energy (x-axis) and of the intensity of 
the VIS dressing field (y-axis); a) experimental results (OD); 
b) theoretical prediction; c) detail of the theoretical prediction 
(notice the extended intensity range). The experimental and 
theoretical results are in good qualitative agreement at low 
intensities. While the experimental signal seems to fade at 
I=4TW/cm^, the calculations predict ongoing resonant fea- 
tures up to intensities as high as 10 TW/cm^. In particular, 
the characteristic flip in the asymmetry of the resonant Fano 
profiles doesn't occur just once, as observed in the experi- 
ment; instead, it is periodic in the intensity of the dressing 
field. See text for more details. 



CT/ ~ 20% /, respectively) is found to smear out the reso- 
nant features at intensities higher than 5 TW/cm^. The 
ac-Stark shift of singly-excited states in helium are known 
to approach the ponderomotive energy [25j , which is the 
limiting value expected for a Rydberg electron bound to 
an unpolarizable core [s^- For the highest terms of the 
autoionizing spj„ series, instead, our calculations show 
that the ac-Stark shift is significantly larger than such 
limit. In fact, in DES the two electrons are strongly cor- 
related: for the higher terms in the 5^2™ series, e.g., the 
parent ion is permanently polarized towards the outer 
electron. This is further indication that traditional for- 
mulas for estimating polarizability, tunneling rates and 
ac-Stark shifts are inapplicable to DES. At higher intensi- 
ties, not shown here, we observe a multi-branch splitting 
of the SP22 doubly-excited state that we associate to the 
multi-photon Rabi oscillation between the spjjj the 2p^ 
^S and the sp2Q states. 

Some simplifying assumptions in the calculations are 
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likely responsible for the occasional discrepancies be- 
tween theory and experiment. In particular, the the- 
oretical XUV pulse has a purely Gaussian spectrum, 
while the experimental spectrum shows some residual 
odd-harmonic modulation. Furthermore, we do not rep- 
resent the wide pedestal that is known to accompany the 
main VIS dressing pulse in practice. Finally, we didn't 
take advantage of the systematic relative uncertainty, of 
the order of ±20%, introduced by the calibration of the 
VIS intensity absolute value. Even so, the remarkable 
agreement with cither the data recorded at the nominal 
intensity of 3.5 TW/cm^ reported in Q and those used 
in the present work is evidence that our conclusions arc 
robust and sound. 

In conclusion, in this work we showed how the inver- 
sion of asymmetric resonant transient absorption profiles 
in ATAS of dressed helium is periodic in the intensity of 
the driving laser, and that the field- free dark ^S state 
actually gives rise to a distinct signal in the ATA spec- 
trum under the infiucncc of a moderately strong VIS laser 
field. The relation of TAS spectra to the non-diagonal 
component of the electrical susceptibility highlighted in 
this work offers an alternative and rational perspective on 
ATAS spectroscopy and suggests a direct way to measure 
and thus modulate the ac-Stark shift in transiently bound 
states, which is an important tool for the implementation 
of quantum-control protocols. Finally, the observation 
that the ac-Stark shift of autoionizing states close to the 
N=2 threshold is larger than the theoretical limit in the 
single-active electron approximation is an indication that 
common model tunncling-ratc formulas arc inapplicable 
to such states and, by extension, to any multi-electron 
atom where polarization of the core affects the effective 
field strength experienced by the ionizing active electron. 
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